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The physical kinetics of single domain ferroelectric ma- 
terials are studied using Landau-Khalatnikov equation. The 
hysteretic curves are obtained numerically. The effective co- 
ercive electric field theoretically varies with the driving ampli- 
tude and frequency. The effects of thermal noise are explored 
using the Fokker-Planck kinetic equation. The ferroelectric 
switching times are discussed and Quantum effects are briefly 
explored. 

PACS: 77.80.Dj, 77.80.Fm, 77.22.Gm 



I. INTRODUCTION 

Ferroelectric materials are endowed with a thermody- 
namic spontaneous electric dipole moment per unit vol- 
ume P s at a minimum of the free energy per unit volume 
F(P, T). Both experimentally |lHlq| and theoretically 



the polarization can be changed in direction via 
a dynamic switching process accompanied by hysteresis 
Hps - 47 1 ■ Some of the partially understood experimen- 
tal properties of ferroelectric domains (27]^-|58| include 
the following: (i) No unique coercive field is evident dur- 



ing the hysteretic process of switching. [E8 -39 (ii) When 
an alternating electric field is applied, a hysteretic loop is 
observed only within a limited frequency band and a lim- 
ited range of the amplitude-frequency product, [fj 43 4(i}| 
(iii) Formation of multiple hysteresis loops (4?]] are com- 
monly observed. 

From the thermodynamic viewpoint, one often as- 
sumes a phenomenological form for the free energy per 
unit volume F(P,T) and employs the thermodynamic 
law 



dF = -SdT + E • dP. 
The thermodynamic electric field is given by 



E 



dF 
&P 



(equilibrium) . 



(1) 



(2) 



The phenomenological free energies adequately describe 
the equilibrium experimental properties of ferroelectric 
materials. 

The dynamical hysteretic properties of the ferroelectric 
domains are less well studied from a theoretical viewpoint 
PHpO -42|. In principle, a changing polarization is equiv- 
alent to a current density J = (dP/dt). One expects a 
current density to give rise to a non-equilibrium electric 
field ^dynamical = = p{dP / dt) where p is an "internal 
resistivity" . The above argument leads to the Landau- 
Khalatnikov dynamical equation of motion |59|-]62| , 



„ fdF\ fdP\ , . , 

\dP ) ^ \~dt ) y non ~ equilibrium). (3) 

The Landau-Khalatnikov Eq.(3) has been employed 
most often in the linearized form close to thermal equi- 
librium 



5E 



■4 up - 



(4) 



where x is the isothermal electric susceptibility. From the 
linearized Landau-Khalatnikov Eq.(4), one may deduce 
the relaxation time 



PX- 



(5) 



While the relaxation time r can be employed to estimate 
the switching time, it is evident that hysteretic loops are 
inherently non-linear. Our purpose is to solve the non- 
linear Landau-Khalatnikov Eq. (3) with a view towards a 
deeper theoretical understanding of ferroelectric hystere- 
sis. The solutions to be discussed in the work which fol- 
lows have been obtained numerically. The qualitative fea- 
tures deduced from the non-linear Landau-Khalatnikov 
Eq.(3) are in satisfactory agreement with experiment. 

In Sec. II, the Landau-Khalatnikov equation is con- 
sidered for the case in which the driving electric field 
has the form E(t) = Eocos(wt). The numerically ob- 
tained hysteretic curves will be exhibited in Sec. III. In 
Sec. IV the effects of thermal noise are explored using 
the Fokker-Planck physical kinetic equation |33| for the 
probability distribution. The ferroelectric switching time 
is computed as a function of the applied electric field. In 
Sec. V we present the Lagrangian form of the equations 
of motion and extend the theory to include second deriva- 
tive terms in the differential equation of motion. The 
quantum mechanical consequences of the Lagrangian for- 
mulation of the physical kinetics are discussed in Sec. VI. 
Future prospects for theoretical extensions are discussed 
in the concluding Sec. VII. 



II. LANDAU-KHALATNIKOV EQUATION 

To measure Hysteretic cycles one applies an electric 
field to the ferroelectric sample E(i) = EoCOs(cL>t). The 
Landau-Khalatnikov equation 



P 



dP 
~dt 



dF(P,T) 
dP 



E cos(wt), 



(6) 



then determines how the polarization P(t) responds in 
time. For a stable hysteretic curve, one looks for a peri- 
odic solution 
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P(t + (27r/w)) = P(i) (cyclic process). 



(7) 



If the polarization P = (0, 0, P) and the applied electric 
field E = (0,0, E) are along an "easy axis", then Eq.(6) 
reads 



(dP\ (dF{P,T) 



M 



dp 



En cos(wi). 



(8) 



We choose for a model free energy, the Landau phe- 
nomenological form 



F(P,T) 



1 



(P 2 -P s {Tf 



(9) 



where P s (T)and Xo{T) are, respectively, the polarization 
and electric susceptibility with a zero static electric field. 

For the purpose of numerical integration of Eqs.(8) and 
(9), one may introduce the dimensionless quantities 



6 = cut, x= (P/P s ), and y = ( X oE/P s ). (10) 



Also 



7] = 



1 



Eqs.(8)-(11) now read 
dx 



and z 



= ( 



En 



\pujP s 



— + rjx(x 2 - 1) = z cos 9. 



(11) 



(12) 



To calculate a hysteretic loop one must numerically solve 
Eq.(12) subject to an initial condition x(8 = 0). After 
integrating Eq.(12) over one period of 2ir, one finds that 



x(9 = 2n)=F(x(9 = 0); 77,20) ■ 



(13) 



The initial condition x{9 = 0) = x for a hysteretic loop 
x{9 + 2tt) = x(9) must be a fixed point of the function 
T; i.e. 



x = T(x-.rj, z ) 



(14) 



Such fixed points exist only for a limited region of values 
of the parameters (77, zq). Thus, only for finite ranges of 
frequency u and driving electric field E in Eq.(8) will 
give rise to clean periodic hysteresis loops. For some 
regimes of frequency and driving electric field, the motion 
in the phase plane (E, P) will be chaotic. The limited 
range of parameters (77, z ) for observing hysteretic peri- 
odic motions in the phase plane is in agreement with well 
known features observed in ferroelectric experiments. 



III. NUMERICAL RESULTS 

Employing Eqs.(l) and (9), we find that the thermo- 
dynamic electric field implied by the model under con- 
sideration, 



(dF\ 
—J (equilibrium), (15) 



obeys 



E = ^7-51 (P 2 p s) (thermal). (16) 



,2X0^ 

The purely thermodynamic Eq.(16) gives rise to hys 
teretic behavior at a switching polarization of 

which occurs at a coercive electric field of 

Ps 



Pswitch 



(17) 



(Pswitch \ 
3xo J ' 



(18) 
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FIG. 1. Shown is a metastable thermal switching hystere- 
sis curve calculated from Eq.(16). If the applied field starts 
negatively at E < —E c and the applied field is subsequently 
quasi-statically increased, then the polarization follows the 
lower curve and jumps to the upper curve when the coercive 
field E c is exceeded. 




(XoE/P s ) 

FIG. 2. The above plot is a numerically simulated hysteric 
loop calculated from the Landau-Khalatnikov model. The 
parameters chosen were rj = 50.0 and zq = 200.0. 



2 



When a cyclic process is numerically simulated for the 
case of finite frequency as in FIG. 2, the resulting hys- 
teretic loop looks qualitatively similar to the thermody- 
namic FIG.l. 

However, there are quantitative features which depend 
on the magnitude of the frequency and the amplitude 
of the applied oscillating electric field. For example the 
dynamical "coercive field" deduced from the numerical 
simulations at finite frequency is quite different from the 
static thermodynamic "coercive field" . Furthermore, the 
"shape" of the hysteretic loop depends strongly on the 
product of the frequency and the internal resistivity via 
2i] = (xopu)^ 1 . This point is illustrated in FIG. 3. 
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FIG. 3. Two different numerically simulated hysteretic 
loops are shown. The smaller (inner) loop corresponds to 
the values r) = 5.0 and zo = 10.0. The larger (outer) loop 
corresponds to the values r\ — 5.0 and zo — 50.0. The two 
loops correspond to the same frequency. However, the larger 
loop is induced by the larger amplitude of the applied electric 
field. 



It is evident, that the dynamically induced coercive 
field (for fixed frequency) depends on the amplitude of 
the applied electric field. Thus, the large (outer) hys- 
teretic loop in FIG. 3 has a much larger dynamic coercive 
field than does the small (inner) hysteretic loop. It is 
not in general possible to deduce the thermodynamic co- 
ercive field of Eq.(18) by the observation of a single dy- 
namical hysteretic loop. From an experimental point of 
view, this ambiguity in determining the thermodynamic 
coercive field is a long standing problem. 



IV. THERMAL NOISE 

The existence of an internal resistivity p in the Landau- 
Khalatnikov Eq.(3) requires a thermally fluctuating elec- 
tric field AE(i); i.e. Eq.(3) with the inclusion of electric 
field noise reads 



dP 

~di 



= E 



dF 
dP 



+ AE. 



(19) 



We have included "white noise" field fluctuations obeying 

'k B T' 



AE(t)AE(tf) 



pl5(t-t'). 



(20) 



If W(P, t)d 3 P denotes the probability for P G d 3 P, then 
the electric field noise in Eqs.(19) and (20) give rise to 
the Fokker-Planck equation 

8W d \(dG\ u .fk B T\dw\ 

p ^r = dp-\{dp) E , T w+ {—)wj> (21) 

where 

G(P,E,T) = P(P,T)-E-P. (22) 
In thermal equilibrium, the probability is determined by 

, N ffl(A(E,T) -G(P,E,T))\ , 

W eq (P) = cxp i V \ bT M . (23) 

where A(E, T) is determined by the total probability nor- 
malization 



W eq (P)d 3 P = 1. 



(24) 



If the polarization P = (0, 0, P) and the applied electric 
field E = (0,0, E) are along an "easy axis", then it is 
sufficient to use a one-dimensional probability W(P, t)dP 
for P G dP which obeys a one-dimensional Fokker-Planck 
equation 



dW d 



'dG 



W 



dt dP } \ dP J ET " \ ' \ n J dP 2 



k B T\ d 2 W 



■ (25) 



The "jump time" associated with Eq.(25) may be com- 
puted in a well known manner. Suppose that the polar- 
ization at time zero is given by P(t = 0) = P,. With a 
final polarization Pf > P,, let r[P/, Pi] be the mean "first 
passage time" for the polarization to move from a region 
P < Pf to the complimentary region P > Pf. In detail, 
let t* be the smallest positive time for which P(t*) = Pf. 
The mean first passage time is then r = t*. 

The closed form expression for the mean first passage 
time is well known to be 

r r A G{p ' E ^- G{p, ' E ^)i kBT dP'dP (26) 

JPi J-oc 

In principle, it is possible to solve the equation 

t = T[P f ,Pi\ (27) 
for the final polarization 
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P f = P(t;Pi). 

Eqs.(26), (27) and (28) imply 
'dP\ ( dF^ 



dt 



dP 



= E + E, 



where 



dP 
~dt 



E = 



k R T\ fd 2 P 



n 



dt 2 



(28) 



(29) 



(30) 



In the limit of large single domains f2 — > oo, we have 
E — > so that Eq.(29) becomes equivalent to the Landau- 
KhalatnikovEq.(3). When \E\ « \E\, then the Landau- 
Khalatnikov equation is sufhciently accurate for comput- 
ing first passage times. However, very near critical tem- 
peratures and or critical electric (coercive) fields, a finite 
domain can switch due to thermal noise. The integral in 
Eq.(26) for the first passage time must then be evaluated. 



V. LAGRANGIAN FORMULATION 

In order to understand the extended Landau- 
Khalatnikov-Tani equation |64| for ferroelectric single do- 
mains, let us first review the extended Landau-Lifshitz- 
Gilbert equations |35|,[36| for ferromagnetic single do- 
mains. For the ferromagnetic case, there exists a free 
energy per unit volume -F(M, T) which obeys the ther- 
modynamic law 

dF = -SdT + H-dM. (31) 

The thermodynamic magnetic intensity is given by 

(32) 



(OF \ 
dM] T ^ uilib ™)- 



For non-equilibrium situations in an applied magnetic 
intensity H, there is an effective field 



H 



eff 



fdF\ fdM\ 



(33) 



The equation of motion for the magnetization then reads 



fdM\ TT 
I I = 7M x H. 



V dt 



l eff, 



(34) 



where 7 is the gyromagnetic ratio. 

For the case of a Ferroelectric, Eqs.(l) and (2) replace 
Eqs.(31) and (32), and the effective magnetic intensity 
Eq.(33) is replaced by the effective electric field 



E e// = E - 



OF 
dP 



dP 

~dt 



(35) 



While in the Landau-Khalatnikov theory, the effective 
field E e // is set to zero, the Tani extension to the 



Landau-Khalatnikov theory includes an inertial term of 
the form 



4n 



d 2 p 



(36) 



Let us consider Eq.(36) in more detail. The dipole mo- 
ment of N charged particles in a volume f2 is defined 

as 



JV 



(37) 



where rj is the position of the j particle with charge 
ey. From the equation of motion for the j th particle with 



mass m 



j ■ 



(38) 



and Eq.(37) one easily obtains Eq.(36) wherein uj p can 
be identified with the plasma frequency 



(39) 



Eqs.(35) and (36) yield the Landau-Khalatnikov- Tani 
theory for ferroelectric domains; i.e. 




4tt\ /d 2 P 



p 



dP\ fdF 



dt \dP 



E. (40) 



Eq.(40) may be written in terms of the effective La- 
grangian L = flC; 



2tt\ fdP 



dt 



C = -3 Kr ) -F(P,T) + E-P (41) 



and the Rayleigh dissipation function R = WZ; 

2n = p {m 



We have 



d ( dC 



- (—\ - ( dU \ 



dt\d(dP/dt)J \dP J \d(dP/dt)J 



(42) 



(43) 



The canonical "momentum" II conjugate to the polar- 
ization P may be defined by 



dC \ _ /4tt^\ fdP 
U= - l \d(dP/dt))-{^ p ){dt 



(44) 



Employing the Lagrangian viewpoint, one may quantize 
the single ferroelectric domain theory in the usual man- 
ner. 
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VI. QUANTUM KINETICS 

In the fully quantum mechanical theory one expects 
the canonical commutation relation 



[n,p] = -mi 



(45) 



which may also be proved employing a microscopic view. 
From Eq.(37) and (44) it follows (with Tj = Vj) that 



n 



4 \ N 



From the microscopic commutation relation 



in 



1, 



(46) 



(47) 



together with Eqs.(37), (39) and (46) follows the macro- 
scopic quantum mechanical commutator in Eq.(45). One 
representation of the quantum mechanical operator II is 
defined by 



n 



-iti 



d_ 



(48) 



Forming the Hamiltonian from the Lagrangian in the 
standard canonical formalism, 

QHe// =n-P-Q£, (49) 

yields the effective Hamiltonian for a ferroelectric domain 



sin 



eff 



Up 



LT 



n {F(P, T) - E • P} . (50) 



where Eqs.(41), (44) and (49) have been invoked. 
Eqs.(48) and (50) describe a single domain ferroelectric 
grain in macroscopic quantum mechanical form. 

From Eqs.(9), (48) and (50) one may compute the 
quasi-classical approximation for the quantum tunneling 
switching time. At zero electric field E — 0, the quantum 
transition rate per unit time to switch polarizations (say 
P s -> -P s ) is given by 



7 exp(-B ) 



quantum 

with a barrier factor given by 




and an attempt frequency of 



7o ( , , 

7T J V "Wo 

In Eq.(53), the Barrier energy per unit volume is 



(51) 



(52) 



(53) 



U Q = 



and 



Thus 



P 2 

± S 

8xo 



4ttxo 



32 /SlUa 



(54) 



(55) 



(56) 



Quantum tunneling through the barrier will dominate 
classical thermal activation over the barrier if 



exp(-Bo) » exp(-QU /k B T). 
Thus for temperatures small on the scale of 



(nuo_\ 

\k B B J 



( 3^wo \ 



(57) 



(58) 



quantum tunneling will dominate classical thermal acti- 
vation. For typical laboratory samples T << T* so that 
the life-time of domains for very small volume domains 
is dominated by quantum mechanics. Nevertheless, for 
large volume domains there is stability and the kinet- 
ics are well described by the Landau-Khalatnikov kinetic 
equation. 



VII. CONCLUSION 

We have shown that physical kinetics, via the Landau- 
Khalatnikov equation, adequately describes laboratory 
hysteretic behavior for large volume single domain fer- 
roelectric materials. The important physical quantities 
include the amplitude and frequency dependence of the 
effective coercive electric field. The formalism describ- 
ing thermal noise as well as quantum noise have been 
explored. Both forms of noise are unimportant for large 
volume ferroelectric domains. In the small domain vol- 
ume limit, quantum noise effects dominate the classical 
thermal noise effects at laboratory temperatures substan- 
tially below the critical temperature. 

The technical importance of ferroelectric switching 
times for small ferroelectric domain volumes is well doc- 
umented |I|-|lS|] . The computation of switching times re- 
quires mesoscopic quantum mechanical models. Thus far, 
the theory of such quantum switches is in its infancy. 
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